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lower  than  1  nni^/Hz  at  1  Hz  and  a  fairly  smooth  decay  prorortictial  to  . 
f~'*  per  decade  up  to  40  Hz,  leading  to  a  value  of  10“^  nirr/Hz  at  10  Hz 
and  10~5  nm^/Hz  at  20  Hz.  Diurnal  variations  due  to  industrial  noise 
sources  such  as  saw  mills  are  mostly  found  in  the  4  Hz  to  10  Hz  band. 
Compared  to  Scandinavia,  the  noise  level  in  Central  Europe  is  ccmpcurable 
or  even  lower  for  frequencies  about  1  Hz  but  it  is  about  a  decade  higher 
at  10  Hz  and  above. 

The  spatial  correlation  properties  of  noise  and  signals  were  investigated 
with  data  from  a  teitporary  9-element  array  with  a  diameter  of  3  km.  In  the 
1  Hz  -  2  Hz  band,  the  noise  coherence  values  drop  to  0.25  at  interstation 
distances  of  800  m,  and  in  the  2  Hz  -  4  Hz  band  at  interstaticn  distancies 
of  400  m,  respectively.  There  is  no  clear  evidence  for  a  negative  noise 
correlation  distance  at  these  frequencies. 

The  relative  amplitudes  of  regional  phases,  i.e.  Pn,Pg,Sn,  and  Lg  are 
strcngly  path-dependent.  From  northern  azimuths  Pq-  and  are 

dominant  with  an  optimum  signal- i-o-noise  ratio  in  the  frequency  band  belcw 
10  Hz.  For  southern  and  south-western  bearings,  vhich  correspcnd  to  Alpine 
wave-paths,  Pn  and  Sn  shew  the  largest  amplitudes  with  a  very  high  frequen¬ 
cy  content  up  to  50  Hz. 

Signal  coherence  in  the  same  frequency  bands  as  they  were  used  for  noise 
correlation,  proved  to  be  excellent  (  >  0.9  )  for  all  regional  phases  over 
the  3  km  aperture  of  the  test-array. 

Ihe  final  configuration  of  the  new  regional  array  in  the  BF-area  -  named 
GERESS  -  has  been  selected  as  a  slightly  enlarged  NORESS-type  ring  geemetry 
with  an  innermost  radius  of  200  m  and  a  largest  radius  of  about  2000  m. 


Accession  Pot 

NTIS  GRA&I  J 
DTIC  TAB  ^ 
Unannounced  [ 
Just  if i } on 

□ 

■ - 4 

By. 


.teSPECTEU 


Table  of  Contents 


page 

I.  Design  and  Siting  of  a  New  Regional 

Seismic  Array  in  Central  Europe  i 

II.  Module  Specifications  and  Hardware  Components  24 

III.  Event  recordings  from  the  Geress  test-array  32 

IV.  Shallow  Refraction  Seismic  Survey  5} 

V.  Acknowledgements  57 


ill 


I.  Design  and  Siting  of  a  New  Regional  Seismic  Array  in  Central  Europe 


1.  Introduction 

Since  tne  establishment  of  the  Graefenberg  (GRF)  array  in  1980,  array  research  in  Ger¬ 
many  was  focused  on  teleseismic  data  due  to  the  broadband  instrument  response  with  a 
high-cut  filter  at  5  Hz  and  an  average  station  distance  of  more  than  10  km.  In  recent  years, 
an  increased  intctcst  emerged  in  the  evaluation  of  high-frequency  signals  from  events  at 
regional  distances.  Small-size  and  many  element  arrays  like  NORESS  have  proven  capable 
of  taking  advantage  of  the  very  efficient  propagation  of  high-frequency  seismic  phases. 
Together  with  a  similar  recently  installed  array  in  Northern  Norway  (ARCESS),  the  pro¬ 
posed  array  in  Southern  Germany  (tentatively  named  GERESS  =  GERman  Experimental 
Seismic  System)  would  built  an  array  network  with  nodes  at  roughly  1000  km  intervals. 
Such  a  configuration  is  particularly  suited  for  automated  location  procedures  applying 
antenna  theory  and  it  is  relevant  in  view  of  a  possible  deployment  of  "in-country"  seismic 
systems  for  a  surveillance  of  a  possible  future  nuclear  test  ban  treaty. 

In  Central  Europe  a  siting  survey  has  carefully  to  investigate  noise  sources  due  to  the 
dense  population,  high  traffic,  and  various  industries.  Apart  from  the  local  noise  conditions, 
sufficient  seismic  signals  from  different  azimuths  should  be  recorded  and  evaluated  to  assess 
the  propagation  of  high-frequency  signals,  especially  across  tectonic  boundaries  like  the 
Tomquist-Teisseyre  suture  between  the  Russian  Platform  and  Western  Europe  or  the 
influence  of  the  Alps  on  signals  originating  from  the  Mediterranean  earthquake  region. 

After  describing  the  site  survey  in  chapter  2,  a  typical  noise  power  spectrum,  taken  from 
data  recorded  in  Southern  Germany,  will  be  compared  with  spectral  noise  conditions  in 
Scandinavia.  A  low-noise  profile  is  only  one  condition  for  a  reasonable  array  site.  Equally 
important  is  the  transfer  function  of  the  receiver  crust  which  influences  the  signal  charac¬ 
teristics  and  the  signal  to  noise  ratio.  Unlike  the  fairly  homogeneous  Scandinavian  shield  the 
geologic  conditions  in  Central  Europe  are  rather  complex  reaching  from  thick  alluvial  sedi¬ 
ments  in  the  North  German  Basin  which  is  separated  by  the  Variscan  mountain  front  from 
the  Alpine  foreland  to  the  Alps  itself.  Besides  the  Black  Forest  at  the  Eastern  Rhinegraben. 
the  Bavarian  Forest  (BF)  area  at  the  border  to  Czechoslovakia  and  Austria  represents  the 
largest  outcropping  crystalline  block  in  Germany.  This  region  is  principally  suited  for  an 
array  installation  of  the  proposed  kind  as  has  been  demonstrated  by  the  excellent  detection 
capabilities  of  conventional  seismic  stations  in  Austria  and  Czechoslovakia. 

In  chapter  3  signal  and  noise  correlation  measurements  in  the  BF-area  will  be  described 
which  are  used  as  important  constraints  for  the  array  design.  For  this  purpose  a  nine-element 
array  was  operated  temporarily  around  the  site  where  the  lowest  noise  values  were  found 
previously. 
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From  these  measurements  the  final  geometry  of  the  GERES  S  array  was  determined  and  its 
array  response  is  compared  with  the  existing  Scandinavian  arrays. 

A  preview  on  the  installation  schedule  of  the  new  array  will  be  given  in  chapter  4  includ¬ 
ing  information  on  the  type  of  instruments  and  the  data  acquisition  units  used. 


2.  Site  Survey. 

Field  work  has  been  carried  out  initially  in  October  1987  and  then  continuously  from 
April  to  July  1988.  Some  additional  data  have  been  collected  later  in  1988  to  prove  the 
long-term  variability  of  noise  conditions  and  to  calibrate  the  results  with  data  from  well- 
defined  events  (e.g.  JVE-explosions  in  Nevada  and  Kasachstan). 

Unlike  NORESS  which  is  built  within  a  subarray  of  the  large  teleseismic  NORSAR  array, 
it  was  not  possible  to  collocate  the  new  high-frequency  array  with  an  existing  GRF-subarray 
because  the  latter  teleseismic  broadband  array  is  placed  on  a  sedimentary  column  which 
strongly  attenuates  high  frequencies.  As  shown  in  Figure  1,  the  new  array  site  is  situated 
about  150  km  east  of  the  GRF-array  location. 

The  noise  measurements  concentrated  on  the  outer  Bavarian  Forest  just  east  of  a  major 
fault  line  named  the  "Pfahl"  which  separates  the  Bohemian  Massif  from  the  western 
Molasse  -  a  quaternary  and  tertiary  sedimentary  basin  in  the  Alpine  foreland.  The  advantage 
of  the  BF-area  is  its  geological  setting  (crystalline  outcropping  rocks)  and  the  low  popula¬ 
tion  density.  The  landscape  is  mostly  mountainous  up  to  1200  m  elevation.  Nearly  all 
recording  sites  were  situated  in  extensive  forest  areas  to  minimize  cultural  noise  and  instru¬ 
ments  were  installed  on  granite  or  gneiss  rocks  to  record  high  frequencies,  especially  from 
events  at  regional  distances.  Due  to  the  installation  on  the  surface,  the  seismometers  were 
quite  sensitive  to  wind  noise  but  for  technical  and  financial  reasons  no  other  arrangement 
was  feasible.  Detection  capabilities  derived  from  noise  estimates  of  these  data  should 
represent  conservative  values. 

For  survey  purposes  three  portacorders  with  direct  recording  were  used.  In  case  of  favour¬ 
able  places  one  of  three  digital  data  acquisition  systems  was  installed  for  a  time  period  of 
several  weeks.  The  digital  systems  were  PCM  recording  instruments  of  Lennartz  5800  type 
with  the  following  specifications: 

ADC  -  66  dB  resolution,  gain  ranging  -  126  dB  dynamic  range, 
sampling  frequency  -  250  Hz,  and  low-cut  filter  -  44  Hz,  6  pole  Bessel. 

Each  PCM  system  was  equipped  with  three  vertical  short-period  seismometers  (1  Hz 
Geotech  S-13)  which  were  installed  at  distances  between  100  m  and  300  m  to  avoid  false 
alarms  by  coincidence  triggering. 

Very  soon  the  site  survey  concentrated  on  an  area  which  is  closest  to  the  CSR  border. 
Very  few  roads,  low  population  density,  and  extensive  woodland  offered  adequate  pre- 
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conditions  for  seismic  installations.  The  first  station  was  established  in  April  1988.  It  was 
placed  on  weathered  granite.  Later,  two  other  stations  were  established  on  gneiss.  During  the 
JVE-experiment  an  additional  station  was  installed  on  granite  at  a  place  where  adequate 
housing  and  power  facilities  were  available  to  operate  a  hub  for  a  temporary  small  array 
subsequently.  All  four  sites  were  located  within  a  radius  of  5  km.  Portacorder  records 
showed  a  generally  low  noise  background  and  especially  a  small  day-to-night  variation.  The 
only  obvious  disadvantage  appeared  to  be  some  saw-mills  which  generated  monochromatic 
seismic  noise  during  working  hours. 

For  comparison  some  noise  samples  have  been  taken  at  the  Graefenberg  array  (station  B5) 
which  confirmed  earlier  measurements  showing  relatively  high  cultural  noise  levels. 

Finally  during  a  short  trip  to  Norway  in  October  1988,  some  recordings  were  made  at 
NOR  ESS  to  get  a  direct  comparison  of  noise  values  by  using  the  same  data  acquisition  sys¬ 
tem  and  -  more  important  -  the  same  processing  procedure  as  for  data  from  the  area  under 
investigation  in  Germany. 

The  processing  of  field  data  included  several  steps.  At  the  beginning  the  PCM  field  tapes 
had  to  be  converted  to  standard  9-track  IBM-compatible  format. 

Having  recorded  the  field  data  in  event  mode  the  pre-event  window  can  be  used  for  noise 
evaluation.  There  are  many  well-established  methods  to  estimate  the  power  spectral  density 
of  stationary  time  series  (e.g.  Oppenheim  and  Schafer,  1975;  Welch,  1967).  In  this  repon 
the  following  procedure  was  applied  ;  Altogether  20  seconds  of  the  pre-event  window  were 
divided  into  19  blocks  of  2  seconds  length  each  with  an  overlap  of  one  second.  Each  data 
block  was  padded  with  zeros  to  get  a  FFT-transformation  length  of  2'®  reminding  that  the 
original  sampling  frequency  of  the  field  data  was  250  Hz.  The  19  raw  Fourier  spectra  were 
averaged  to  lower  the  variance  without  affecting  stationary  noise  peaks.  The  final  step 
includes  an  average  of  12  noise  spectra  from  each  station  and  a  plot  of  the  mean  values  and 
their  standard  deviation. 

The  noise  spectra  were  calculated  separately  for  day  and  night.  They  cover  the  whole  time 
period  during  which  the  corresponding  station  was  operating.  By  that  procedure,  working 
hours,  weekends,  and  different  weather  conditions  are  included  in  the  noise  estimate. 

Figure  2  shows  a  typical  noise  spectrum  for  the  BF-area.  The  solid  line  represents  the 
average  and  the  dotted  lines  one  standard  deviation  calcuiated  as  explained  in  the  preceding 
text.  The  estimate  includes  data  from  summer  and  from  winter  when  the  ground  was 
covered  with  snow.  In  the  individual  spectra  no  significant  difference  could  be  observed. 
The  spectrum  shows  a  continuous  slope  with  a  small  variance  and  the  2  Hz  peak,  commonly 
observed  at  stations  in  Central  Europe,  can  only  be  recognized  during  night  times.  The  day 
spectmm  is  dominated  by  a  noise  maximum  at  4-5  Hz  which  is  supposed  to  originate  from 
a  saw  mill  at  a  distance  of  a  few  kilometers.  Apart  from  this  peak  the  noise  spectrum  shows 
a  smooth  decay  proportional  to  f  from  1  Hz  to  30  Hz  and  a  small  standard  deviatir 
Taking  into  account  that  the  recording  time  covered  nearly  a  period  of  four  months,  this 
area  seemed  to  be  promising  for  an  array  installation. 
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To  secure  this  suggestion,  three  other  sites  were  explored  in  the  vicinity.  Indeed,  their 
spectra  looked  very  similar  and  confirmed  the  favorable  opinion  about  the  area. 

To  put  our  noise  results  from  the  Bavarian  Forest  into  proper  place,  additional  data  were 
collected  at  the  Graefenberg  array  and  at  NORESS.  GRF  can  be  regarded  as  a  typical 
Western  European  site  whereas  NORESS  is  well  known  for  its  excellent  noise  conditions 
being  situated  on  the  Scandinavian  shield. 

At  day  the  GRF-spectrum  (Figure  3)  shows  much  larger  noise  values  than  in  any  of  the 
BF-spectra.  The  influence  of  industry  and  traffic  is  especially  pronounced  as  seen  in  the 
large  variation.  The  night-spectrum  at  GRF  looks  similar  to  the  BF-spectra  between  1  Hz 
and  8  Hz.  For  higher  frequencies  the  different  geological  setting  (sediments)  causes  still 
higher  noise  values.  From  this  direct  comparison  we  can  conclude  that  the  BF-area  exceeds 
most  other  places  in  the  TRG  -  and  certainly  GRF  -  as  a  potential  site  for  an  establishment 
of  a  high-frequency  array. 

More  interesting  is  the  comparison  of  the  proposed  Bavarian  area  with  NORESS.  The 
spectra  shown  in  Figure  4  were  calculated  from  a  24  hour  noise  sample  analyzed  with  the 
same  procedure  as  described  above.  There  are  remarkable  differences  in  the  noise  spectra. 
For  low  frequencies  around  1  Hz  NORESS  clearly  suffers  from  the  influence  of  the 
Norwegian  coast  which  results  in  an  order  of  magnitude  higher  PSD-values  compared  to  the 
Bavarian  Forest  area  (lOnnAHz  to  less  than  Inm^/Hz)  .  Apparantly  these  microseisms 
lead  to  a  steep  slope  of  the  spectrum  proportional  to  f  up  to  frequencies  of  2  -  5  Hz. 
For  higher  frequencies  this  slope  is  flattened  and  becomes  comparable  to  the  f  fall-off  at 
the  BF-area.  The  absolute  noise  values  at  10  Hz  and  20  Hz  are  certainly  lower  at  NORESS 
but  also  some  influence  of  industrial  noise  between  5  Hz  and  8  Hz  can  be  identified  in  the 
spectrum.  According  to  personal  information  from  the  NORSAR  staff,  there  is  also  a  saw 
mill  operating  at  a  distance  of  roughly  10  km  from  the  NORESS  array. 

In  summarizing  the  comparison  between  NORESS  and  the  BF-area  we  found  higher  noise 
values  at  NORESS  for  frequencies  below  2  Hz  and  higher  noise  values  at  the  BF-site  for 
frequencies  between  2  Hz  and  20  Hz.  Tne  consequence  of  this  difference  in  terms  of  detec¬ 
tion  capabilities  can  only  be  evaluated  comparing  events  recorded  at  both  arrays. 


3.  Array  Design. 

The  GERESS  array  is  planned  as  part  of  a  multi-array  network  which  will  include  the 
NORESS  and  ARCESS  arrays  in  Scandinavia.  Consequently  the  design  philosophy  fol¬ 
lowed  those  existing  arrays  and  it  is  well  documented  in  several  publications  (e.g.  Followill 
and  Harris,  1983,  Mykkeltveit,  1985), 
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The  geometry  of  these  regional  arrays  is  based  on  concentric  rings  spaced  at  log-periodic 
intervals  in  radius  R,  following  the  relation 


R  =  a",n  =  0,  1.2,  3 


(1) 


with  3,  5,  7,  and  9  elements  in  each  ring,  plus  one  in  the  center. 

With  an  odd  number  of  elements  in  each  ring,  the  corresponding  coarray  samples  the  local 
wavefield  in  the  least  redundant  way.  In  this  sense,  such  a  configuration  comprises  subsets 
of  sensors  with  very  different  intersensor  separations  implying  that  both  high-frequency  and 
low-frequency  regional  phases  can  be  well  enhanced  by  appropriate  subsets  of  the  array 
(Kvaerna,  1989). 

An  essential  aspect  of  the  array  design  is  the  supposition  that  intersensor  spacings  should 
correspond  to  maximum  correlation  for  the  signal  and  minimum  correlation  for  the  noise. 
Seismic  signal  and  noise  characteristics  are  influenced  not  only  by  the  geologic  formations 
of  the  array  site,  but  also  by  the  type  and  strength  of  sources  and  their  locations  with  respect 
to  the  array.  Consequently,  determination  of  the  actual  values  of  the  parameters  in  (1)  for  a 
specific  site  is  to  be  based  on  observed  signal  and  noise  correlations. 

Our  correlation  curves  are  derived  from  measurements  made  with  a  9-sensor  test-array 
operating  from  December  1988  to  April  1989  in  the  BF-area  where  we  found  the  lowest 
noise  conditions  earlier.  The  configuration  of  the  test-array  and  its  corresponding  coarray  are 
shown  in  Figure  5.  Interstation  distances  range  from  200  m  to  3000  m.  Short-period  instru¬ 
ments  (Geotech  S-13)  were  placed  in  small  vaults  of  2  m  to  4  m  depth.  At  this  depth  the 
overburden  of  soil  was  penetrated  in  general,  but  the  underlying  crystalline  bedrock  was  still 
heavily  weathered.  The  main  advantage  of  the  vault  installation  was  to  avoid  wind  noise 
induced  by  the  surrounding  forest. 

New  digital  data  acquisition  systems  (Geotech  PDAS-100)  were  used  for  recording  which 
sampled  at  100  Hz  with  a  resolution  of  16  bits.  Data  synchronization  was  controlled  by  an 
external  radio  signal  (DCF-77).  The  tcsi-array  operated  in  detector  mode  with  a  conventional 
STA/LTA-detector,  recording  length  was  generally  3  min  with  a  pre-event  window  of  30 
sec. 

Figure  6  shows  an  average  noise  power  spectrum  of  the  test-array.  The  trend  and  absolute 
values  coincide  quite  well  with  the  results  of  the  previous  site  survey  (Figure  2)  confirming 
confidence  that  the  selected  area  is  suited  for  an  array  installation.  Noise  correlations  were 
calculated  from  the  pre-event  window  of  35  different  events  covering  the  whole  recording 
period.  The  raw  data  were  bandpass  filtered  between  1  Hz  -  2  Hz  and  2  Hz  -  4  Hz,  respec¬ 
tively  with  a  two-pole,  both-directional  (zero-phase)  Butterworth  filter. 


The  correlations  are  expressed  by 


X;  (  x(i)  -  X  )  (  y(i)  -  y  ) 

1  E  (  x(i)  -  5f  I  [  X  (  y(i)  -  y  1 

i  =  1  i  =  1 

x(i)  and  y(i)  are  sample  values  for  sensor  x  and  y,  N  is  the  number  of  samples  and  x  and  y 
the  mean  values  of  the  N  samples. 

The  window  length  was  taken  as  10  sec  (1000  samples).  The  correlations  were  calculated 
for  all  36  different  combinations  of  seismometer  pairs,  i.e.  all  the  coarray  points  (Figure  5). 

Our  results  are  summarized  in  Figure  7.  In  the  lower  Jf,  the  individual  correlation  values 
are  shown  within  the  1  Hz  -  2  Hz  passband  together  with  a  smootned  average.  The  large 
scatter  in  the  data  may  reflect  as  well  the  short  window  length  as  also  some  problems  with 
the  synchronization  during  data  acquisition.  On  the  other  hand,  the  solid  line  representing  an 
average  trend,  indicates  a  quite  reasonable  behaviour,  especially  when  we  add  the  upper  half 
of  the  figure  which  presents  t.he  corresponding  data  in  the  2  Hz  -  4  Hz  passband.  Defining 
the  correlation  length  as  that  station  separation  for  which  the  noise  coherence  drops  below 
0.25,  we  get  a  correlation  length  of  about  800  m  for  the  1  Hz  -  2  Hz  passband  and  400  m 
for  the  2  Hz  -  4  Hz  passband.  In  a  wide  sense,  this  result  is  consistent  with  values  found  for 
NORESS.  The  main  difference  to  noise  correlation  studies  in  Scandinavia  is  the  absence  of 
a  significant  minimum.  The  noise  correlation  curves  for  the  BF-area  show  a  smooth  decay 
without  a  pronounced  negative  correlation  distance.  If  this  preliminary  result  will  be 
confirmed  by  a  comprehensive  analysis  of  the  future  array  data  there  could  be  various  expla¬ 
nations;  firstly,  the  isotropic  noise  model  (Backus  et  al,  1964)  might  not  be  applicable  in  an 
azimuthally  heterogeneous  geologic  region  like  Central  Europe  or  alternatively,  the 
minimum  in  noise  correlation  for  Scandinavian  stations  originates  from  directional  propagat¬ 
ing  noise  which  could  be  higher  mode  microseisms.  The  higher  noise  level  at  low  frequen¬ 
cies  which  was  earlier  mentioned  in  comparing  Figures  2  and  4  supports  the  latter 
hypothesis. 

However,  in  the  final  design  of  the  NORESS  array,  not  much  use  has  been  made  of  the 
noise  minimum  because  in  the  range  of  interesting  frequencies  the  corresponding  interstation 
distance  varies  and  optimum  geometries  would  be  vastly  different  (Mykkeltveit,  1985). 

As  a  consequence  of  the  noise  correlation  curves,  the  minimum  station  distances  for  the 
planned  GERESS  array  could  be  slightly  increased  in  comparison  to  NORESS  taking  into 
account  that  the  detection  window  with  the  best  signal  to  noise  ratio  is  generally  shifted  to 
lower  frequencies  for  stations  in  Central  Europe.  Whether  this  can  lead  to  a  larger  overall 
aperture  depends  on  signal  correlation.  Figure  8  shows  a  regional  earthquake  recorded  by 
the  test-array  in  the  BF-area.  It  is  a  mining-induced  event  originating  at  Lubin  in  the  Polish 
copper  mines  at  a  distance  of  360  km  from  the  test-array.  Besides  the  first  arriving  Pn- 
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phase,  we  recognize  Pg  and  Lg  wavetrains  as  dominant  signals,  a  typical  picture  for  events 
from  northern  or  nonhwestem  azimuths. 

In  Figure  9  the  signal  spectra  are  plotted  together  with  the  background  noise.  Although  all 
three  phases  stay  above  the  noise  level  up  to  frequencies  of  30  Hz,  the  best  signal  to  noise 
ratio  is  obviously  in  the  1  Hz  -  5  Hz  band.  Certainly  this  figure  changes  with  azimuth  and  at 
the  BF-site  there  are  special  wave  paths  along  the  Alps  where  we  see  much  higher  signal 
frequencies  and  Pn  and  Sn  become  the  dominant  phases.  But  in  general,  frequencies  below 
10  Hz  are  most  important  for  the  evaluation  of  regional  waves  within  1000  km  distance. 

Signal  correlations  were  calculated  with  the  same  procedure  which  was  used  to  establish 
the  noise  correlation  curves.  Time  shifts  have  been  performed  with  an  optimum  line-up  of 
peaks  and  troughs  at  the  different  sensors,  and  1  sec  or  2  sec  windows  of  each  phase  were 
analyzed.  In  Figure  10  the  signal  correlation  curves  are  shown  for  the  Pn,  Pg,  and  Lg  phases 
using  the  same  1  Hz  -  2  Hz  and  2  Hz  -  4  Hz  passbands  as  for  the  noise  correlation.  Addi¬ 
tionally  the  noise  correlation  is  marked  in  Figure  10  by  stars  which  represent  the  plus/minus 
one  standard  deviation  range.  As  can  be  seen  from  the  figure,  all  signal  correlation  values 
stay  above  0.9  for  all  interstation  distances  out  to  3000  m  with  the  exception  of  the  Lg 
phase  in  the  2  Hz  -  4  Hz  passband.  Similar  signal  coaelation  curves  were  obtained  analyz¬ 
ing  a  number  of  regional  events  from  different  azimuths. 

FTom  these  measurements  it  was  concluded  to  maintain  the  general  geometry  of  the 
NORESS  array  but  to  enlarge  it  for  the  planned  GERESS  array  by  a  factor  of  4/3  which 
results  in  a  radius  of  200  m  for  the  innemtost  ring  and  about  2000  m  for  the  outermost  ring, 
respectively.  In  conclusion,  the  parameters  in  equation  (1)  for  the  concentric  array  layout 
were  specified  as  =  200  m  and  a  =  2.15  . 

The  final  siting  of  the  25  array  elements  had  to  take  some  local  geologic  and  topographic 
peculiarities  into  account.  In  addition  some  restrictions,  imposed  by  the  forest  authorities 
had  to  be  considered.  Nevenheless  the  final  configuration,  displayed  in  Figure  11  is  still 
sufficiently  close  to  the  concentric  ring  concept  and  consequently  the  GERESS  array  will 
represent  a  uniform  element  in  the  multi-array  network  including  ARCESS  and  NORESS, 
the  arrangement  of  the  latter  is  displayed  in  Figure  12  (left  side). 

For  comparison  the  array  responses  of  GERESS  and  NORESS  are  plotted  on  the  righ 
side  of  Figure  11  and  Figure  12,  respectively. 

Due  to  the  larger  aperture,  the  main  lobe  is  smaller  than  the  NORESS  main  lobe  but  on 
the  other  hand  minor  side  lobes  can  be  seen  in  the  reject  region  in  the  wavenumber  domain. 
For  a  fixed  number  of  sensors  this  reflects  the  wellknown  trade-off  between  resolution  and 
aliasing  propenies.  Whether  the  superior  resolution  of  the  GERESS  array  will  lead  to 
improved  location  and  phase  identification  capabilities  has  to  be  po.stponed  until  the  new 
array  becomes  operational.  An  essential  aspect  of  the  array  design  using  the  concentric  ring 
concept  deals  with  the  uniform  average  of  the  coarray  domain.  As  shown  in  Figure  13,  the 
distribution  of  coarray  points  is  very  similar  for  GERESS  and  NORESS  and  different 
subgeometrics  will  achieve  the  optimum  signal  to  noise  gain  in  different  frequency  bands 


and  with  different  wavelengths. 
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4.  Discussion  and  Preview. 

Since  several  conflicting  demands  can  be  put  on  array  performance,  it  is  difficult  to  make 
general  comparisons  of  array  design.  For  the  regional  array  concept  developed  at  NORSAR, 
the  underlying  objectives  are  determined  in  the  context  of  a  monitoring  system  for  a  future 
nuclear  test  ban  treaty.  A  special  aspect  of  such  a  system  is  the  identification  of  weak 
seismic  events  using  high-frequency  data  (Archambeau  et  al,  1986).  As  a  prototype  array 
NORESS  has  demonstrated  promising  capabilities.  It  is  however  difficult  to  generalize 
NORESS  results  because  this  array  is  located  within  on  old  .shield  region  with  high  Q 
values.  Additionally  the  natural  seismicity  within  regional  distances  is  relatively  weak.  The 
GERESS  array,  on  the  other  hand,  will  be  situated  in  a  rather  complex  geologic  and  geotec- 
tonic  environment  where  high  Q  and  low  Q  wave  paths  drastically  change  with  azimuth. 
The  array  site  is  surrounded  by  various  artificial  seismic  sources  like  quarries  and  mines  but 
more  important,  the  Alpine  earthquake  belt  and  the  Mediterranean  earthquake  zone  lie 
within  about  1000  km  distance.  Consequently,  the  number  and  nature  of  regional  seismic 
events  recorded  at  GERESS  should  be  more  manifold  than  at  any  recording  site  in  Scandi¬ 
navia.  For  optimum  performance  of  a  seismic  surveillance  system,  detection,  location,  and 
identification  of  events  at  low  signal  to  noise  ratios  is  important.  It  will  be  interesting  to 
observe  whether  automatic  procedures  developed  for  these  tasks  at  NORSAR,  can  be 
adapted  to  the  seismic  situation  in  Central  Europe.  In  any  case,  the  use  of  arrays  as  com¬ 
ponents  of  a  global  seismic  network  can  be  evaluated  more  realistically  the  more  the  arrays 
and  their  recordings  represent  typical  conditions  of  a  global  monitoring  system. 

In  summer  1990  the  GERESS  array  should  become  fully  operational.  After  the  final 
configuration  was  defined,  excavation  of  vaults  and  trenching  of  cables  started  in  October 
1989.  Figure  14  shows  that  in  addition  to  the  25  vertical  instruments  four  sites  will  include 
horizontal  components.  All  the.se  are  short-peritxl  (1  Hz)  Geotech  GS-13  type  instruments 
sampled  at  40  Hz.  The  array  will  be  supplemented  by  one  three -component  set  of  GS-13’s 
sampled  at  120  Hz  and  these  instruments  will  be  collocated  with  a  three-component  set  of 
broadband  seismometers  (BB-13)  sampled  at  10  Hz.  The  sites  were  selected  to  have  an 
option  for  a  good  three-component  subarray  taking  into  account  the  site-specific  noise  condi¬ 
tions.  In  contrast  to  vault  A2,  the  central  vault  NU  did  not  reach  the  basement  rock  and  it 
appeares  to  be  noisier  than  the  average  array  station.  Site  C2  is  very  close  to  the  hub  facility 
and  allows  easy  access  which  might  be  favorable  in  case  of  modifications  at  the  broadband 
or  high-frequency  instruments.  As  all  sites  arc  equipped  with  data  acquisition  units  using  24 
bit  A/D  converters,  the  data  will  provide  sufficient  resolution  and  dynamic  range  to  evaluate 
small-size  and  large-size  earthquakes  simultaneously.  For  interesting  events,  data  from  the 
three-component  high-frequency  and  broadband  instruments  can  be  pieced  together  to  study 
the  seismic  wavefield  in  a  very  broadband  sense. 
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'ITie  GERESS  array  will  be  installed  and  operated  as  a  joint  project  of  the  Geophysical 
Laboratory  of  Southern  Methodist  University  Dallas  and  the  Geophysical  Institute  of  Ruhr- 
University  Bochum.  Data  will  be  available  at  NORSAR  and  Bochum. 
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Geographical  location  of  the  existing  teleseismic  GRF-broadband 
array  and  the  planned  regional  GERESS  array. 
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Fig.  2 


:  Displacement  power  spectra  of  noise  in  the  Bavarian  Forest  (BF)  area  for  day 
and  night.  Solid  line  represents  an  average  of  12  different  data  samples, 
dotted  lines  mean  plus/minus  one  standard  deviation 
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Fig.  3  :  Displacement  power  spectra  of  noise  at  the  GRF-array  (site  B5). 

For  explanation  see  Fig.  2. 
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Fig.  4  :  Displacement  power  spectra  of  noise  at  NORESS  (site  C4), 
For  explanation  sec  Fig.  2. 
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Fig.  5  :  Nine-element  test-array  in  the  BF-area.  Coordinates  of  central  vault:  48.845N,  13.700E 
The  coarray  on  the  right  side  shows  all  interstation  distances  in  vector  space. 
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F’.g.  6  :  Average  displacement  power  spectra  of  noise  for  all  elements  of  the  test-array 
in  the  BF-area.  Dotted  lines  mean  plus/minns  one  standard  deviation. 


correlaUoQ 


-  15  - 


noise  correlelion.  f  =  2  —  4  IIz 


Fig.  7  :  Noise  corrclarion  values  within  the  test-array,  (top  :  prefilter  2  Hz  -  4  Hz, 
bottom  :  I  Hz  -  2  Hz  passband)  The  solid  lines  represent  a  smoothed  average  curves. 
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Fig.  8  :  Test-array  data  (9  traces)  of  an  ml=3.8  (KBA)  earthquake  near  Lubin  (Poland). 
Epicentral  distance  is  354  km.  Prominent  regional  phases  arc  marked  as  Pn,  Pg,  and  Lg. 


Fig.  9  :  Signal  and  noise  power  spectra  for  Lubin  event  for  site  A 1  from  Fig.  8. 
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e084.3:  signal/noise  correlations,  f  =  1  Hz  —  2Hz 


intersensor  distance  [m] 


e0a4.3:  signal/noise  correlations,  fsSHz  —  4  Hz 


intersensor  distance  [m] 

Poland  25.  03.  1989  A  =  354  km  mL  =  3.8 

Fig.  10  :  Signal  coherence  curves  (top:  1  Hz  -  2  Hz  passband,  bottom  :  2  Hz  -  4  Hz  passband) 
for  the  test-array  calculated  from  data  of  the  Lubin  event  shown  in  Fig.  8.  The  stars 
represent  plus/minus  one  standard  deviation  of  the  preceding  noise  window. 
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Fig.  1 1  :  Final  configuration  of  the  GERESS  array  including  25  vertical  seismometer  sites. 
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Fig.  13  :  Comparison  of  NORESS  (lop)  and  GERESS  (bottom)  coarray  points 
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Fig.  14  ;  The  geometry  of  the  GERESS  array  including  station  names  and  instrumentation 
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II.  Module  Specifications  and  Hardware  Components 
developed  for  the  GERESS  Test-Array. 


1 .  Functional  Requirements 

To  get  an  reliable  estimate  of  the  capability  of  the  planned  array,  the  noise  and  signal 
coherence  had  to  be  determined  over  short  term  and  seasonal  changes.  This  resulted  in  a 
projected  measurement  time  of  about  half  a  year  and  a  measurement  area  covering  the 
whole  aperture  of  the  planned  array. 

The  calculation  of  signal  and  noise  coherency  in  dependence  of  station  distance  requires  a 
long  term  timing  accuracy  of  1  msec  between  all  stations. 

So  the  overall  design  requirements  for  the  test  array  installation  were  : 

9  stations  distributed  over  the  array  aperture  resulting  in  station  distances  between 
300m  -  2200  m  from  the  central  station 

timing  accuracy  of  1  msec  in  the  whole  network 
sampling  rate  of  100  Hz 

digital  data  acquisition  with  high  dynamic  range  and  resolution 
measurement  time  of  half  a  year  -  mostly  in  winter  1988/89 

unattended  operation  of  the  remote  stations  for  at  least  one  week  because*  of  the  lim¬ 
ited  personnel 

event  oriented  data  storage  because  of  the  huge  amount  of  raw  data  produced  in  high 
frequency  mode 

2.  Available  equipment 

At  date  of  the  planned  test  installation,  the  only  digital  field  measurement  units  available 
off  shelf  were  Geotech’s  PDAS  100,  in  it’s  first  production  run  with  rudimentary  software, 
however.  They  fulfilled  the  requirements  of  high  precision  digital  data  acquisition,  but  failed 
in  some  other  important  requirements  listed  above: 

Timing  accuracy  of  free  running  units  was  excellent  on  a  short  term  base,  but  unac¬ 
ceptably  bad  for  a  long  term  installation. 

Power  consumption  would  require  a  change  of  the  battery  pack  every  few  days. 

The  limited  memory  and  lack  of  a  coincidence  detector  to  avoid  false  alarms  would 
even  in  an  event  triggered  mode  require  data  transfer  to  a  laptop  on-site  every  few 
hours. 
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I'he  test  array  was  to  be  installed  in  a  forest  area  with  few  roads,  the  operators  were 
expected  to  be  quite  busy  with  processing  and  archiving  data  in  the  central  station  at  the 
SUN'  workstation.  So  any  solution  requiring  frequent  on-site  inspections  was  impractical. 
Instead  power  supply,  timing  synchronization,  coincidence  trigger  and  data  transfer  from 
PDAS  to  PC  had  to  be  perfonrted  from  the  central  station.  Fortunately  the  array  site  allowed 
for  long-line  cable  installation  and  the  electronic  laboratory  of  the  Institute  of  Geophysics 
had  to  build  the  necessary  interface  modules. 


3.  System  Configuration. 

The  overall  technical  design  of  the  test  array  is  shown  in  Fig.  15  : 

Eight  remote  units,  consisting  of  one  S13  and  PDAS  100  each,  are  connected  via  4-pair 
telephone  cable  with  the  central  station.  One  pair  is  for  power  supply  (by  60V  DC  and 
DC/DC  conversion  to  12V),  another  for  1  PPS  timing  signals.  The  remaining  pairs  form  a 
short  haul  modem  connection  between  PDAS  and  PC  (RD  and  TD  only,  software 
handshake),  while  the  hardware  handshake  (DCD,  DTR)  is  used  for  trigger  in/out. 

The  central  unit  was  connected  to  a  three-component  set  of  S-13  seismometers.  It  evaluates 
the  trigger-in  from  all  stations  by  its  digital  10  facility  and  provides  the  coincidence  trigger 
for  the  whole  array  (trigger-out).  Additionally  a  DCF-receiver  provides  the  1  PPS  signals  for 
the  whole  array.  Data  transfer  is  done  by  connection  of  laptop  2  with  one  remote  PDAS  via 
serial  link  (KERMIT).  During  this  time  (.5  MB  in  approx.  20  min)  data  acquisition  has  to 
be  suspended  due  to  the  single  task  architecture  of  the  current  software  version  in  PDAS. 

Figure  16  is  the  circuit  diagram  of  the  driver  card  in  the  central  station  rack  for  each 
remote  unit.  The  key-switch  disconnects  the  power  supply  and  initiates  a  remote  boot-up  of 
the  PDAS.  The  BC  177  is  the  driver  for  the  current  loop  in  the  1  PPS  transmission. 

Figure  17  shows  the  pulse  generator  for  the  1  PPS  by  a  monofiop  74121  and  amplifiers 
for  each  remote  unit. 

Figure  18  gives  details  of  the  connection  from  PDAS  unit  0  (the  central  unit)  and  laptop  2 
to  the  remote  units  via  modems  1-8  (type  MA-I4/MA-18  of  Hedin  Tex  AB,  Sweden). 
Selection  of  one  remote  unit  for  data  communications  will  disable  trigger  transmission  via 
IC  4001  for  this  station,  so  the  selection  for  unattended  operation  is  ’O’. 

Figure  19  shows  the  connection  box  at  each  remote  station.  Its  purpose  is  to  provide  the  sig¬ 
nals  needed  by  PDAS  1-8  from  the  connection  to  the  central  station  via  telephone  line. 
Power  is  supplied  by  a  DC/DC  converter,  while  opto-coupler  CNY  17-4  decodes  the  1  PPS 
signal.  Hardware  handshaking  of  the  modem  MA-19  is  used  for  trigger  in/out,  while 
1C  4528  and  relays  K1  and  k2  are  needed  for  the  remote  boot  up  by  disconnection  of  power 
in  the  central  station  (because  of  the  buffered  power  supply  by  the  battery  pack  this  would 
not  work  otherwise). 
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4.  Experiences 

The  half-year  installation  time  proved  most  of  the  above  design  criteria.  The  on-site 
inspections  were  reduced  to  monthly  intervals  needed  only  for  check  of  the  mechanical 
installations.  Data  transmission  and  power  supply  worked  without  difficulty. 

Some  problems  however  arose  due  to  the  immature  software  in  the  PDAS.  Without  any 
systematics  the  clock  sometimes  confused  the  1  PPS  resulting  in  a  1  second  time  offset. 
This  was  seldom  however  and  quite  easy  to  verify  and  to  correct  in  the  data  subsequently. 
Another  problem  was  the  unaccurate  timing  of  the  detector.  Due  to  the  time  windows  pro¬ 
cessed  on  block  by  the  signal  processor  of  PDAS,  detection  time  was  only  a  rough  estimate 
and  dependent  on  the  relative  timing  of  processing  blocks  to  seismic  signal  in  each  station. 
This  demanded  an  unreasonably  long  time  window  for  the  coincidence  detector  not  corre¬ 
lated  to  signal  propagation  time  in  the  array  any  more. 


M.  Joswig 
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Fig.  15  :  Block  diagram  of  GERESS  test  array 


Fig.  17  :  Circuit  dia^',’  m  of  switch  unit  PDAS  1-8  to  laptop  2 


Fig.  18  :  Circuit  diagrani  of  line  driver  for  1  PPS  timing  signal 
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III.  Event  recordings  from  the  Geress  test-array 


As  outlined  in  chapter  II,  the  data  acquisition  system  principally  operated  in  detection 
mode.  A  conventional  STA/LTA  -  defector  was  implemented  on  the  single  PDAS- 100  chan¬ 
nels  and  a  time  window  was  set  in  which  a  minimum  number  of  channels  had  to  trigger. 

Using  an  updating  LTA  the  recording  length  was  variable  but  generally  complete  seismo¬ 
grams  for  regional  events  from  Pn-waves  to  Lg-waves  were  recorded. 

After  in,spection  of  the  detection  times  which  were  compared  with  bulletin  data,  i.e.  PDE 
for  tclcseismic  events  and  EMSC  for  regional  seismicity,  the  corresponding  wavefomis  were 
stored  in  a  data  base  using  CSS  2.8  format. 

The  operation  of  the  test-array  covered  a  time  period  from  December  8  1988  to  April  4 
1989.  Due  to  technical  problems  related  to  the  temporary  installation,  the  actual  recording 
time  was  only  92  days  of  the  total  time  period. 

Altogether  276  events  were  recorded  and  summarized  in  an  event  list.  112  events  were 
classified  as  local  within  an  epicenter  distance  of  150  km,  they  originated  mostly  from 
quarry  blasts  in  Czechoslovakia  and  Germany.  Additionally  67  teleseismic  events  (  distance 
greater  than  2000  km  )  were  recorded  and  finally  97  events  occurred  at  regional  distances  ( 
150  km  -  2000  km  ). 

As  the  planned  array  is  specifically  aiming  at  high-frequency  signals  particular  emphasis 
was  given  to  the  evaluation  of  these  regional  events.  Table  1  contains  a  list  of  nine  events 
which  were  selected  for  even  azimuthal  coverage.  The  epicenters  of  these  events  are  shown 
in  Figure  20.  We  start  our  discussion  with  event  1  (  Figure  21  )  originating  from  Lubin  in 
Poland  which  was  already  displayed  in  Figure  8.  In  a  previous  report  (  Haijes,  1989  )  it  has 
been  mentioned  that  two  mining  districts  in  Poland  cause  remarkable  seismicity,  namely  the 
copper  mine  near  Lubin  and  the  coal  mining  district  near  Katowicz. 

These  events  exhibit  very  consistent  Pg-  and  Lg-phascs  which  dominate  the  seismograms. 
A  clear  Pn-onset  is  often  missing  in  the  seismograms  from  the  coal  mines  and  this 
difference  can  be  attributed  to  the  different  geology  in  tl.e  source  region. 

Seismic  waves  from  event  2  (  Figure  22  )  traveled  the  shortest  distance  of  265  km  from 
the  epicenter  in  Austria  to  the  GERESS  test-array.  With  its  sharp  Pn-phase  the  event  looks 
explosion-like  although  it  is  reported  as  a  shallow  earthquake  and  a  close  look,  indeed, 
shows  a  weak  starting  phase.  About  4  sec  after  the  first  arrival,  we  see  a  clear  P*-phase 
whereas  the  Pg-phase  is  almost  hidden  in  the  coda.  Accordingly,  the  S-wavetrain  starts  with 
the  S*-phase  and  no  distinct  Sg-  or  Lg-wave  is  visible. 

Continuing  further  to  the  South  on  our  epicenter  map,  event  3  (  Figure  23  )  shows  an 
earthquake  from  Hungary  with  an  emergent  Pn-phase  followed  by  a  long  coda.  No  clear 
later  arrivals  can  be  recognized  although  the  Lg-amplitudes  are  easily  detected. 
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The  following  two  events  (  no  4  and  no  5  from  table  1  )  which  are  plotted  in  Figure  24 
and  Figure  25  originate  from  the  same  source  region  at  the  Adriatic  coast  of  Jugoslavia. 
Ti'.ey  set  an  e.xample  for  the  differences  in  phase  amplitudes  which  can  exist  for  very  close 
epiccnters.  While  event  4  contains  a  clear  Pg-phase  but  almost  no  Lg-phase,  event  5  shows 
just  the  opposite.  In  addition,  the  Sn-phase  is  also  very  different  for  both  earthquakes>  event 
4  exhibits  two  distinct  S-phases,  the  first  of  which  is  certainly  the  Sn-phase,  while  the 
second  could  be  interpreted  as  an  S*-phase.  Consequently,  these  two  earthquakes  would 
constitute  a  challenging  task  for  an  automated  event  processor. 

The  next  three  events  (nos  6  -  8  ),  plotted  in  Figures  26  -  28,  come  from  southeastern 
azimuths  and  the  seismic  waves  travel  Alpine  wave  paths.  The  amount  of  high  frequencies 
in  all  three  recordings  is  surprising,  Pn-  and  Sn-  waves  are  dominant,  both  pha.ses  are  fol 
lowed  by  a  long  coda.  Unfortunately  the  recording  window  for  two  of  these  events  was  too 
short  to  discuss  the  Lg-wave  amplitudes.  For  event  6,  originating  in  Northern  Italy,  the  data 
quality  is  poor  and  only  three  traces  are  recorded  for  the  expected  arrival  time  for  Lg-waves 
but  there  appears  to  be  no  signal.  The  situation  is  similar  to  the  earlier  discussed  events 
from  Jugoslavia  where  Lg-waves  appear  and  disappear  within  small  epicenter  differences. 
The  same  observation  holds  for  the  Pg-wave.  From  a  clear  onset  in  event  7  (  Figure  27  ) 
the  Pg-amplitude  decreases  to  the  Pn-coda  level  for  event  8  (  Figure  28  ). 

The  last  regional  earthquake  from  table  1  is  plotted  in  Figure  29  (  event  9  ).  Waves  from 
this  event  strike  the  test-array  from  a  northw'estern  azimuth  and  the  general  picture  of  the 
seismogram  looks  like  event  1  which  had  a  bearing  to  the  northeast.  The  crustal  phases  Pg 
and  Lg  dominate  the  seismogram,  the  Pn-phase  shows  very  small  amplitudes  and  Sn  is  not 
di.scemible. 

In  concluding  the  description  of  regional  seismograms  recorded  at  the  test-array,  two  limi¬ 
tations  should  be  mentioned;  firstly,  the  source  mechanism  and  size  of  the  events  will  also 
influence  the  relative  excitation  of  different  wave  types  secondly  th**  ^'^'aracterization  of 
regional  seismograms  by  the  few  traditional  phases  can  only  be  regarded  as  a  first  order 
approximation. 

Nevertheless,  the  differences  in  regional  phases  can  also  be  demonstrated  in  the  frequency 
domain.  In  Figures  30  -  32  a  power  spectrum  of  those  phases  together  with  the  background 
noise  is  plotted.  Figure  30  shows  the  main  difference  between  an  event  from  northern 
azimuth  (  no  1,  Poland  )  and  an  earthquake  from  southwestern  azimuth  (  no  8,  Switzerland). 
In  both  cases  the  sampling  rate  was  100  Hz,  i.e.  the  antialiasing  filter  cuts  the  information  at 
40  Hz.  Whereas  the  spectral  amplitudes  of  the  signal  'bases  for  the  Polish  event  stay  above 
the  noise  spectrum  for  frequencies  below  20  Hz  with  an  optimum  SNR  around  2  -  3  Hz,  the 
Swiss  earthquake  shows  excellent  SNR’s  out  to  40  Hz. 

It  would  certainly  have  been  interesting  to  record  even  higher  frequencies.  This  unusual 
frequency  behaviour  is  confirmed  by  the  second  Swiss  earthquake  (  no  7  in  Figure  32 
which  strikes  the  test-array  from  roughly  the  same  bearing.  Due  to  the  lower  sampling  rate 
of  50  Hz  the  antialiasing  filter  acts  already  at  20  Hz. 
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Events  from  southeastern  azimuths  (  Figure  31,  events  no  3  and  no  4  )  show  a  rather  steep 
spectral  slope  comparable  to  the  Polish  event  but  in  contrast  to  that  event  the  Pn  -  phase  of 
the  Hungarian  and  Yugoslavian  earthquake  exhibit  the  highest  frequencies. 

Regarding  the  high  frequencies  which  are  a  common  feature  of  regional  events  at  distan- 
cies  below  1000  km,  a  sampling  rate  of  40  Hz  -  as  planned  for  the  final  array  design  -  is  to 
be  seen  as  a  minimum  requirement.  Higher  sampling  rates  would  be  desirable  and  the  one 
high-frequency  three-component  seismometer  which  is  planned  to  sample  at  120  Hz  will 
become  of  great  importance. 

At  distances  of  1000  km  and  more  the  waveforms  from  earthquakes  in  Greece  and  Turkey 
as  well  as  from  southern  Italy  and  northern  Africa  are  very  complex,  no  significant  signal 
energy  can  be  recognized  above  10  Hz.  This  ob.servation  is  in  accordance  with  earlier  results 
from  the  site  survey. 

For  the  same  set  of  regional  events  the  location  capabilities  of  the  test-array  were  investi¬ 
gated.  As  has  been  shown  earlier  (  compare  Figure  10  ),  the  signal  correlation  proved  to  be 
excellent  over  the  aperture  of  the  test-array.  On  the  other  hand,  a  conventional  narrow-band 
wavenumber  analysis  resulted  in  large  variances  for  the  azimuth  i,stimates  depending  on  the 
passband  of  the  frequency  prefilter.  Figure  33  shows  the  k-analysis  for  the  Polish  event  (  no 
1  )  for  four  different  frequencies.  The  estimated  bearing  fluctuates  from  34  degrees  at  1  Hz 
to  18  degrees  at  4  Hz.  In  this  situation  a  new  method  of  direction  estimation  of  wideband 
signals  (  Nawab  et  al,  1985  )  gave  excellent  results.  The  method  uses  the  zero-delay  covari¬ 
ances  of  the  complex  analytic  representation  of  the  seismometer  outputs  to  estimate  the 
wavenumber  spectrum.  The  power  along  each  radial  direction  in  this  wavenumber  spectrum 
is  proportional  to  the  temporal  power  spectrum  of  plane  wave  sources  in  that  direction.  Fig¬ 
ure  34  gives  the  wavenumber  plots  for  the  Pn  phase  cf  6  events  from  table  1.  The  estimated 
bearings  match  the  theoretical  values  given  in  table  1  quite  well.  If  one  takes  into  account 
that  the  geometry  of  the  test-array  was  far  from  optimal,  the  location  capabilities  of  the  final 
GERESS  array  should  allow  an  independent  evaluation  of  the  seismicity  in  Central  Europe. 
Locations  derived  from  GERESS  data  can  be  used  as  starting  solutions  in  a  location  pro¬ 
cedure  including  single  stations  from  a  European  network. 
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region 

date 

origin  time 

lat. 

Ion. 

A/km 

az. 

1.  Poland 

25.03.1989 

12:46:40.3 

51.62°N 

16.irE 

354 

28° 

2.  Ausnna 

11.02.1989 

02:46:12.2 

47.97°N 

17.03°E 

265 

110° 

3.  Hungary 

27.01.1989 

03:55:05.6 

47.04“N 

16.95°E 

315 

128° 

4.  Yugoslavia 

16.12.1988 

11:35:52.3 

44.72°N 

15.04°E 

470 

167° 

5.  Yugoslavia 

13.02.1989 

00:35:13.0 

44.79°N 

14.69°E 

457 

170° 

6.  Northern  Italy 

25.12.1988 

18:27:36.1 

44.97°N 

9.05°E 

558 

221° 

7.  Switzerland 

07.01.1989 

02:29:42.4 

46.36°N 

7.49°E 

543 

242° 

8.  Switzerland 

02.04.1989 

06:59:00.6 

47.19°N 

9.01°E 

390 

244° 

9.  GDR 

13.03.1989 

13:02:18.2 

50.74°N 

9.84°E 

349 

309° 

Table  1  :  List  of  regional  events  recorded 
with  GERESS  test-array 


DISTANCE  [km] 
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DISTANCE  [km] 


GERESS  TEST-ARRAY 


•  EPICENTER 


Fig.  20  :  Epicenter  map  of  regional  events  recorded  with  GERESS  test-array 


event:  e084.3  dole:  MAR  25  (084),  1989  slarltime:  12:47:08,739 


Fig.  21  :  f  vnt  i  ;  Poland  03/25/1989  12:46:40.3 


event:  e027.l  dole:  JAN  2/  (027).  1989  slorllime:  0J:55:18.199 


727/1989 


event;  e351.1  dole:  DEC  16  (351).  1988  slortlime:  11:36:27.000 


40 


Event  4  :  Yugoslawia  12/16/1988  11:35:52.3 


evenl:  eOlll  dole;  fCO  li  {QH),  1089  slotlllrne:  ()0;35J1.120 


VfRI 


evenl:  eJ60.l  dole:  DEC  25  (360).  1988  slorltime:  18;28:20.'I79 


g.  26  :  Event  6  :  Northern  Italy  12/25/1988  18:27:36.1 


event;  e007.1  dole;  JAN  07  (007),  1989  slofllime;  02:30:^6.790 


43 


Switzerland  04/02/1989 


event:  e092.3  dole:  APR  02  (092),  1989  slorltime:  06:59:80.750 


Switzerland  01/07/1989  02:29:42.4 


event:  e072.1  dole:  MAR  13  (072),  1989  slorllirne:  13:02:42.969 


Fig.  29  :  Ev^nt  9  :  GDR  03/13/1989  13:02:18.2 
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Displacement  power  spectra  of  events  no  1  and  no  8  from  Table  1 
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Fig.  31  : 


Displacement  power  spectra  of  events  no  3  and  no  4  from  Table  1 


Displacement  power  spectra  of  events  no  6  and  no  7  from  Table  1 
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IV.  Shallow  Refraction  Seismic  Survey 


A  detailed  geologic  map  of  the  array  area  is  displayed  in  Figure  35.  A  comparison  with 
topography  shows  that  mainly  the  mountain  tops  consist  of  outcropping  crystalline  rock,  i.e. 
granite  and  gneiss,  whereas  in  the  valleys  a  soil  overburdens  the  basement  rock.  Due  to  this 
heterogeneity,  the  final  siting  of  the  array  vaults  was  very  difficult. 

On  the  one  hand  the  concentric  ring  geometry  had  to  be  bom  in  mind,  on  the  other  hand 
as  many  seismometer  sites  as  possible  should  be  placed  on  bedrock. 

An  extensive  refraction  seismic  survey  was  carried  out  to  measure  the  depth  of  the 
bedrock.  Altogether  43  reversed  refraction  profiles  were  shot  which  were  distributed  over 
the  whole  array  area. 

A  twelve-channel  Geometries  data  acquisition  system  and  a  12  kg  hammer  source  was 
used.  The  spread  of  the  in-line  profiles  varied  between  24  m  and  36  m  (  geophone  distance 
of  2  m  or  3  m  respectively  )  with  an  offset  of  2-3  m  or  12-13  m.  To  increase  the  signal-to- 
noise  ratio,  10-15  hammer  blows  were  stacked. 

The  results  of  the  survey  are  summarized  in  Table  2.  The  station  names  correspond  to 
Figure  14.  Often  more  than  one  profile  had  to  be  measured  at  each  site  and  although  the 
difference  in  location  between  profiles  A.B  or  A,B,C  was  only  of  the  order  of  100  m,  a 
large  variation  in  interval  velocities  was  found.  This  result  reflects  the  complex  surface  geol¬ 
ogy  of  the  area  with  varying  overburden  and  weathering.  In  general  the  refraction  survey 
yielded  a  three  layer  model:  The  first  layer  with  interval  velocities  below  1000  m/s  consists 
of  soil,  the  second  layer  with  interval  velocities  between  1000  m/s  and  2000  m/s  can  be 
interpreted  as  more  or  less  compact  granitic  or  gneissic  mud.  Finally,  the  third  layer  with 
seismic  velocities  clearly  above  2000  m/s  represents  weathered  crystalline  rock.  In  very  few 
cases,  crystalline  bedrock  with  velocities  above  4000  m/s  was  reached  with  the  small-scale 
refraction  method. 

As  a  consequence  of  these  results,  it  is  not  possible  to  place  all  seismometers  of  the  25- 
element  array  on  bedrock,  given  a  target  depth  of  4  -  5  m  for  the  vaults.  If  the  concentric - 
ring  configuration  is  to  be  followed,  the  only  alternative  would  be  the  installation  of 
borehole  seismometers.  As  no  such  instrumentation  is  planned  for  the  GERESS  array,  it  was 
made  allowance  of  a  deviation  of  the  strict  array  geometry.  By  this  procedure,  most  of  the 
array  sites  are  assumed  to  show  noise  values  around  10  nm/s  whereas  a  few  mud  sites  may 
yield  higher  noise  values  up  to  50  nm/s.  Because  these  high  values  originate  from  a  site 
amplification  effect,  also  the  seismic  signals  will  be  increased,  and  therefore  the  signal-to- 
noise  ratio  is  comparable  to  the  bedrock  sites. 

Array  processing  should  take  care  of  these  different  site  characteristics  and  no  ma  i  'c 
disadvantage  in  array  performance  is  expected. 


mica  gneiss 
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mud  (Jliesserde) 
consolidated  mud  (Jliesserde) 


•  GERESS  array  station 


Fig.  35:  Geological  map  of  ihe  GERESS  array  area 
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Table  2:  Results  of  shallow  refraction  seismic 

survey 

Station  name 

Interval  velocities 

Layer  thicknesses 
upshooting  downshooting 

Null  A 

vl»  330  m/s 
v2‘»  1550  m/s 

3.2  m 

3.1  m 

Null  B 

vl-  360  m/s 
v2-  1200  m/s 

4.5  m 

1.7  m 

Null  C 

vl-  370  m/s 
v2-  810  m/s 

1.7  m 

2.2  m 

A3  A 

vl-  380  m/s 
v2»  900  m/s 
v3-  2400  m/s 

2.0  m 

10.0  m 

2.0  m 

8.0  m 

A3  B 

vl-  400  m/s 
v2-  1050  m/s 

3.0  m 

1.3  m 

B1 

vl-  400  m/s 
v2-  900  ra/s 
v3-  2800  m/s 

1.0  m 

12.0  m 

3.0  m 

7.4  m 

B3 

vl-  370  m/s 
v2-  640  m/s 
v3-  1900  m/s 

2.3  m 

9.3  m 

3.0  m 

7.5  m 

B4  A 

vl-  340  m/s 
v2-  950  m/s 
v3-  2000  m/s 

1.4  m 

10.0  m 

2.2  m 

9.3  m 

B4  B 

vl-  400  m/s 
v2-  800  m/s 
v3-  3100  m/s 

1.6m 

4.0  m 

0.5  m 

11.0  m 

0 

B5  A 

vl-  300  m/s 
v2-  850  m/s 
v3-  3000  m/s 

1.0  m 

4.0  m 

0 

B5  B 

vl-  410  m/s 
v2-  1000  m/s 
v3-  2200  ra/s 

1.5  m 

4.0  m 

1.5m 

4.4  m 
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Cl  A 

vl-  370  m/s 
v2-  1400  m/s 
v3-  2500  m/s 

4.0m 

5.0  m 

2.0m 
7.8  m 

Cl  B 

vl-  360  m/s 
v2-  830  m/s 
v3-  2800  m/s 

2.0  m 

7.0  m 

2.0  m 
5.6m 

Cl  C 

vl-  350  m/s 
v2-  700  m/s 

2.0  m 

3.3  m 

Cl  D 

vl-  400  m/s 
v2-  800  m/s 
v3=  1900  m/s 

1.3  m 

6.4  m 

1.0  m 
6.0  m 

C3 

vl-  360  m/s 
v2-  1050  m/s 
v3-  3300  m/s 

2.7  m 

11.0  m 

C4  A 

vl-  600  m/s 
v2-  900  m/s 
v3-  4400  m/s 

2.5  m 
9.4  m 

C4  B 

vl-  420  m/s 
v2-  670  m/s 
v3-  2050  m/s 

1.4  m 

7.0  m 

1.7  m 
1.9  m 

C4  C 

vl-  730  m/s 
v2-  1900  m/s 

5.5  m 

6.8  m 

C5  A 

vl-  370  ra/s 
v2-  1400  m/s 
v3-  3400  m/s 

3.4  m 

9.0  m 

2.3  m 
7.0  m 

C5  B 

vl-  330  m/s 
v2»  900  m/s 
v3=  1600  m/s 

0.5  m 

5.6  m 

1.3m 
2.0  m 

C5  C 

vl-  340  m/s 
v2-  2100  m/s 

2.8  m 

2.4  m 

C6  A 

vl-  370  m/s 
v2-  2500  m/s 

2.0  m 

2.0  m 

C6  B 

vl-  410  m/s 
v2-  1100  m/s 
v3-  2000  m/s 

1.2  m 

5.5  m 

1.0  m 
5.0  m 
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C7 

vl- 

420 

m/s 

3.0 

m 

4.0 

m 

v2- 

2800 

m/s 

10.0 

m 

10.0 

m 

v3- 

5500 

m/s 

D2 

A 

vl- 

360 

m/s 

3.5 

m 

4.7 

m 

v2- 

1900 

m/s 

D2 

B 

vl- 

360 

m/s 

3.0 

m 

5.0 

m 

v2“ 

2100 

m/s 

D2 

C 

vl- 

360 

m/s 

4.7 

m 

2.7 

m 

v2" 

2500 

m/s 

D4 

A 

vl- 

450 

m/s 

1.9 

m 

2.1 

m 

v2- 

1300 

m/s 

7.9 

m 

5.5 

m 

v3“ 

4200 

m/s 

D4 

B 

vl- 

500 

m/  s 

2.2 

m 

11.8 

m 

v2- 

1050 

m/s 

19.6 

m 

7.5 

m 

v3- 

4800 

m/s 

D4 

C 

vl- 

340 

m/s 

1.8 

m 

1.4 

m 

v2- 

1300 

m/s 

D5 

A 

vl» 

360 

m/s 

1.1 

m 

1.2 

m 

v2- 

1000 

m/s 

5.2 

m 

5.0 

m 

v3“ 

1800 

m/s 

D5 

B 

vl- 

340 

m/s 

1.0 

m 

0.1 

m 

v2- 

900 

m/s 

2.5 

m 

1.6 

m 

v3- 

1400 

m/s 

D5 

C 

vl- 

370 

ro/s 

1.4 

m 

1.0 

m 

v2- 

1030 

m/s 

5.6 

m 

3.0 

m 

v3- 

2200 

m/ s 

D6 

A 

vl- 

520 

m/s 

10.1 

m 

2.7 

m 

v2- 

2200 

m/s 

D6 

B 

vl- 

390 

m/s 

1.0 

m 

0.3 

m 

v2- 

1320 

m/s 

15.0 

m 

9.3 

m 

v3- 

5600 

m/s 

D7 

A 

vl- 

350 

m/s 

1.5 

m 

v2- 

600 

m/s 

4.0 

m 

v3- 

4000 

m/s 

D7  B 
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D7  B 

vl-  370  m/s 
v2-  1600  m/s 
v3-  3700  m/s 

1.3  m 
6.3m 

1.5  m 
9.4  m 

D8  A 

vl-  660  m/s 
v2-  1700  m/s 

5.0m 

5.0  m 

D8  B 

vl-  530  m/s 
v2-  1050  m/s 
v3-  2800  m/s 

4.0m 
11.0  ra 

2.0  m 
10.0  m 

D9  A 

vl-  340  m/s 
v2“  700  m/s 
v3-  3400  m/s 

0.7  m 

16.0  m 

1.2  m 
13.0  m 

D9  B 

vl-  340  m/s 
v2“  2200  m/s 

4.3m 

5.5m 

D9  C 

vl-  350  m/s 
v2-  800  m/s 

1.6  m 

3.8  m 
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